Theoretical evidence for strong correlations and incoherent metallic state in FeSe 
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The role of electronic Coulomb correlations in iron-based superconductors is an important open 
question. We provide theoretical evidence for strong correlation effects in FeSe, based on dynamical 
mean field calculations. Our ab initio spectral properties first demonstrate the existence of a lower 
Hubbard band. Moreover, together with significant orbital-dependent mass enhancements, we find 
that the normal state is a bad metal over an extended temperature range, implying a non-Fermi 
liquid due to formation of local moments. Predictions for angle-resolved photoemission spectroscopy 
are made. 

PACS numbers: 71.27.+a, 74.70.Xa, 74.25.Jb 



I. INTRODUCTION 

The discovery of high-temperature superconductivity 
in iron-based compounds triggered a tremendous amount 
of experimental and theoretical research. In general 
terms one can divide these materials into two classes. 
The first is based on iron-pnictogen FePn -1 planes, which 
have to be intercalated with some spacer atoms. The sec- 
ond class is built up by charge-neutral iron-chalcogenide 
FeCn layers, where Cn can be S, Se, or Te. Supercon- 
ductivity in this '11' family was first reported by Hsu et 
al.^ with a transition temperature of T c w 8K for a- 
FeSeo.85- The properties of this material can be further 
modified by pressure^— excess Fe or Sc deficiencies. An- 
other interesting route is alloying FeSe and FeTe yielding 
the family of FeSe^Tei-z superconductors 5 -. 

In contrast to the high-temperature superconduct- 
ing cuprates, whose parent compounds are Mott insu- 
lators, the question about the strength of correlations in 
the iron-based superconductors is not settled. FcSei_ x 
has been studied with angular-integrated photoemission 
(PES) in Refs. 00, while angular-resolved photoemis- 
sion studies (ARPES) have been performed on Fei + . T Tc 8 
and FeSeajTci-z (with x = 0.3 in 0] and x = 0.42 in 
flCjj). For all these compounds, these experiments re- 
veal a significant iron bandwidth narrowing, by a fac- 
tor in excess of ~ 2. Interestingly, recent ARPES ex- 
periments on FeSe^Tei-a; with x = 0.42^ report large 
mass enhancements deduced from the low-energy Fermi 
velocities, ranging from 6 to 23, while a smaller enhance- 
ment (in the range 2 — 3) was reported for Fei+^Te^. 
Specific heat measurements 1 ^ for FeTeo.67Seo.33 gave 
7 = 39mJ/molK 2 , whereas for FeSeo.88 a smaller value 
7 = 9.17mJ/molK 2 was reported* 1 - 

Band structure calculations based on density- 
functional theory (DFT) of FeSe and FeTe have showni 2 - 
that the one-electron band structure of these materials 
is similar to the other iron-based superconductors, as re- 
gards the Fe-d states around Fermi level. However, a re- 



cent first-principles calculation of the screened Coulomb 
interactions gave significantly higher values for the 11 
compounds. 13 Dynamical mean-field (DMFT) calcula- 
tions at T = - albeit within the iterative perturbation 
theory - could reproduce some aspects of the experimen- 
tal results, by taking the interaction as a parameter. 14 
Given the experimental controversies, further insight 
from ab initio theoretical predictions, in comparison with 
experiments, is needed in order to characterize the role 
of correlations for the 11 family relative to other iron 
superconductors (e.g. the 1111 family as LaFeAsO). 

In this article, we investigate the correlation effects 
on the electronic structure of ev-FeSe in the framework 
of DMFT, using numerically exact Monte Carlo simu- 
lations and taking into account the full rotationally in- 
variant Coulomb interactions, evaluated from first prin- 
ciples. We find that this material displays clear evidence 
of strong correlations. First, our results demonstrate 
that the satellite feature observed in PES at a bind- 
ing energy of about —2 eV— >£ should be interpreted as 
a lower Hubbard band of iron origin. Second, we find 
that three of the iron orbitals are characterized by a 
rather low quasiparticle weight, and short quasiparticle 
lifetimes. These features indicate deviation from conven- 
tional Fermi liquid behavior possibly associated with the 
formation of local moments controlled by the Hund cou- 
pling. The possible occurence of local moments in multi- 
band systems was discussed in a model calculation 1 ^ 
and for 1111 pnictide o 16 ' 17 for parametrized interaction 
strengths. The situation in the 1111 is subtle. However, 
taking the interaction strenght calculated from first prici- 
ples, LaFeAsO shows moderate correlations^ Here, on 
the contrary, our present calculations with ab initio inter- 
actions strongly suggest that correlation effects arc most 
clearly revealed in the 11 family. We make theoretical 
predictions for the ARPES spectrum of FeSe. Signifi- 
cant deviations from the DFT band-structure are found, 
which do not simply amount to an overall bandwidth nar- 
rowing. A downward shift of the hole-like bands near the 
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T-poiiit and an upward shift of the electron-like bands 
near the M-point are found. Although ARPES has only 
been reported at this stage for FeSe^Tei-a; , this is con- 
sistent with the trends reported for these alloys.^^ 



II. METHODS 

Our calculations use the recent implementation of 
the combined DFT-DMFT method in a full-potential 
augmented plane wave electronic-structure framework^ 
based on the Wicn2k package.— Localized Wannicr-likc 
orbitals are constructed from an energy window com- 
prising the Fe-3cZ bands and the Se-4p bands, calcu- 
lated in the local-density approximation (LDA). We use 
the tetragonal crystal structure, space group P4/nmm, 
as reported by Margadonna et aZ.j22 using the experi- 
mental value for the Se position in the unit cell. A 
many-body self-energy, computed from DMFT using a 
strong-coupling continuous-time quantum Monte Carlo 
algorithm^ is applied to the subspace spanned by the 
Fe-3ci orbitals. The matrix of Coulomb interaction pa- 
rameters is calculated from first principles using the con- 
strained random-phase approximation (cRPA)^ as ap- 
plied to the iron-based superconductors in Refs. THHl]- 



All screening transitions arc included, except the ones 
within the Fe-3d manifold. In the notation of Refs. [l3ll23j 
this corresponds to the so-called d-dp construction. In 
order to avoid an orbital dependent double-counting cor- 
rection, we do not use the calculated interaction matrices 
directly. Instead, we follow the procedure described in 
Ref. Iia and use the orbitally-averagcd Coulomb inter- 
actions. For FeSe, this yields a local Coulomb integral 
U = F = 4.06 eV and a Hund's coupling J = 0.91 eV. 
This is to be compared with the significantly smaller val- 
ues U ~ 2.7 eV, J ~ 0.8 eV for LaFeAsOJ^^ Please 
not that the definition of J is used here in terms Slater 
integrals, and hence differs from the definition of J in 
Refs. T3II23I I . We will show that the inclusion of the full 



rotationally-invariant Hund's coupling (including spin- 
flip and pair-hopping) is crucial for FeSe. If not otherwise 
specified, all calculations were performed at a temper- 
ature T = 290 K using the fully-localized-limit double 
counting^ Spectra were obtained from the imaginary- 
frequency Monte Carlo data using the stochastic Maxi- 
mum Entropy method.— 



III. RESULTS 

In Fig. [I] we compare the total momentum-integrated 
spectral function obtained within LDA and within 
DMFT. A substantial bandwidth reduction is obtained 
for the Fe-d states near the Fermi energy. Furthermore, 
the Se-p ligand states are shifted to larger binding ener- 
gies as compared to LDA. In addition, there is a structure 
appearing at a binding energy between —1 and — 2eV. In 
order to identify the physical nature of this peak, we 




FIG. 1: (Color online) Comparison of the total LDA DOS 
(black) to the spectral function obtained with DMFT (red). 
'LHB' denotes the lower Hubbard band. The inset shows 
the evolution of the LHB as a function of U, with J scaled 
accordingly (J = 0.9, 1.1, 1.3, respectively.) 
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FIG. 2: (Color online) Orbital-resolved comparison of the 
density of states between the LDA (thin black) and DMFT 
(thick, red) results. The orbitals in the bottom row show 



performed calculations for larger interactions. For com- 
putational efficiency, the calculations for varied U and 
J were done here with density-density interactions only. 
Non density-density terms affect mostly the quasiparticle 
states and little the LHB. It is obvious from the inset of 
Fig.[T]that the peak shifts towards larger binding energies 
with increasing U, which is the expected behavior for a 
lower Hubbard band (LHB), and distinguishes it from a 
low-energy quasi-particle excitation. The existence of a 
lower Hubbard band is nicely consistent with experimen- 
tal PES results. Independent measurements on FeSe&£ 
revealed a broad feature at a binding energy of about 
about — 2eV. It was shown in Ref. Q that the photon- 
energy dependence of this peak (labeled 'B' in Fig. 1 of 
Ref. [a]) indicates that it is of Fe-d origin. Given our 
theoretical results, we propose that the feature observed 
in Refs. @0] in this range of binding energies should be 
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interpreted as a lower Hubbard band. To our knowledge, 
this is indeed the first observation of a lower Hubbard 
band in iron-based pnictide and chalcogenide supercon- 
ductors. Indeed, the shoulder at about -1.5 eV observed 
in PES for LaFcAsO in Ref. [26| can be explained from 
band theory. We also note that the position of the Se-p 
bands agrees well with the experimental peak positions 
of —6 and — 4eV, corresponding to peaks 'C and 'D' in 
Fig. 1 of Ref. @. 

Our results reveal a marked orbital dependence of the 
correlation effects, with significantly stronger correlations 
for the three orbitals d xy and d xzyz , while the d z 2 and 
d x 2_ y 2 orbitals display weaker correlations. Our conven- 
tions for the (x,y) axis are rotated by 45 degrees as 
compared to the crystallographic axis, so that the d xy 
orbital in our definition is the one pointing from Fe to 
Se. The three orbitals displaying stronger correlations 
are thus the ones which form the peak of the LDA-DOS 
(and hence have higher weight at the Fermi level), while 
the two other ones display a 'pseudo-gap' in the LDA 
DOS. In Fig. [2J we plot the spectral functions of each 
Fe-d orbital, calculated for the same parameters as be- 
fore. It is apparent that the LHB discussed above is vis- 
ible mostly in the d xy and d xz ^ yz partial DOS. We have 
also calculated Z m = (1 — J|^) _1 |i W -j-0j which - in a 
Fermi liquid (sec however below) - can be interpreted as 
the quasiparticle weight, and the inverse single-particle 
lifetime Im£ m (zO + ) of each orbital. We find that the 
d z 2 and d x 2_ y 2 orbitals yield larger quasiparticle weights 
Z z 2 = 0.38, Z x 2_ y 2 = 0.47 and longer lifetimes (corre- 
sponding to -Iml](i0 + ) of order 0.04cV at T = 290K). 
In contrast, a linear fit to the slope of the self-energy for 
dxy.xz,yz orbitals would yield lower values of the weights 
Z xy = 0.20, and Z xz yz = 0.28. However, the rather 
short lifetimes (corresponding to — lmY,(i0 + ) xy ~ 0.19 eV 
and —lmTj{iQ + ) xz / yz ~ 0.08 eV) shed serious doubt on 
the validity of a coherent quasiparticle interpretation for 
those orbitals. A further decrease of the temperature 
from T = 290K to T = 190K does not reduce E(i0+). 
In order to identify coherent quasi particles the tempera- 
ture has to be lower than the width of the quasi-particle, 
i.e. T < ZlmS(i0 + ). Here ZImE(i0+) for the d xy or- 
bital corresponds to 440 K, which is higher than 290 or 
190 K and, hence, no coherent quasi-particle can be ex- 
pected. We checked that changing the double counting 
to "around-mean-field"— even increases the incoherence 
(— ImS(i0 + ) X2 .j, z ps 0.4). We elaborate in more detail 
on this unconventional metallic state further at the end 
of this paper. The stronger degree of correlations in 
FeSe, compared to LaFeAsO^, manifests itself also in 
the partial charges of the Fe-d electrons. Due to elec- 
tronic correlations, the electron charge diminishes from 
6.37 in LDA to 6.07 within LDA+DMFT, which should 
be compared to the values for LaFeAsO of 6.40 (LDA) 
and 6.28 (LDA+DMFT). Please note that these changes 
do not correspond to any doping effect, since the total 
electron count of the crystal is still integer. It is merely 
a redistribution of charges due to changes in the hybridi- 
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FIG. 3: Comparison of the momentum-resolved spectral func- 
tion of DFT (top) to the DFT+DMFT results (bottom). 




FIG. 4: Comparison of the low-energy spectral function, pro- 
jected to orbital character, of DFT (left) to the DFT+DMFT 
results (right). Top row: d xy orbital. Bottom row: Degener- 
ate d xz ,d vz orbitals. 



sation between the Fe and the ligand atoms. 

We have also calculated the momentum-resolved spec- 
tral function along high-symmetry directions, displayed 
in Fig. [3] and compared to the LDA bandstructure. 
These results can be viewed as predictions for future 
ARPES experiments on pure FeSe. They can also be 
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compared qualitatively with available ARPES data^ii . 
bearing in mind however that those actually concern 
the Te-substituted compound FeSe^Tei-a;. The overall 
renormalization of the bands and the band- width reduc- 
tion are apparent. The inverse lifetime (ImE) increases 
rapidly with increasing frequency, which results in rather 
broad structures below -0.3 eV. The DMFT quasiparticle 
bandstructure cannot be deduced from a mere rescaling 
of the LDA bandstructure. The strong orbital and en- 
ergy dependence of the self-energy renormalizes the effec- 
tive crystal-field splitting between the orbitals, resulting 
in differential shifts of the quasiparticle bands near the 
Fermi level. An interesting effect is seen for example 
near the T point. In LDA, one can see essentially only 
two dispersing hole-like bands, the outermost band be- 
ing quasi-degenerate. By introducing correlations, this 
degeneracy is lifted, and we can see three hole-like exci- 
tations around the T point. This is qualitatively con- 
sistent with ARPESJ^ We also observe that the hole 
pockets near the T-point are pushed downwards in en- 
ergy by correlation effects, while the electron-like ones 
near the M-point are pushed upwards, in agreement with 
the ARPES spectra presented in Refs. jafiot The or- 
bital character of the LDA+DMFT bands, see Fig. H 
reveals that the outermost band around T and the elec- 
tron pocket around M are dominantly of d xy character, 
which displays the largest effective mass in our calcula- 
tions, see above. This is again in qualitative agreement 
with ARPES, where the heaviest masses were associated 
to the outermost hole pocket and the electron pocket. 
Furthermore, the broader feature at around -0.3 eV is of 
d z 2 character, consistent with experiment. 

Within the DMFT approach used in the present work, 
the effective mass enhancement of electrons in a given 
orbital m is related to the inverse of the weight Z rn calcu- 
lated above. This suggests effective mass enhancements 
in the range ~ 2 for the x 2 — y 2 orbital to ~ 5 for the xy 
one. Direct comparison of these values to experiments is 
difficult in the absence of ARPES data for FeSc. ARPES 
measurements on FeSe^Tei-^^ have been interpreted 
as yielding very large effective masses (m*/mb an d be- 
tween 6 and 23)^°-. This is qualitatively consistent with 
the much larger specific heat coefficient reported for these 
alloys jii in comparison to FeSei It should also be kept 
in mind that the precise extraction of effective masses 
(renormalized low-energy Fermi velocities) from ARPES 
data by comparison to the DFT band structure is sub- 
ject to rather large uncertainties, especially in view of the 
orbital-dependent shifts of the bands induced by correla- 
tions (see above). 

Considering the full rotationally-invariant Hund's 
couplin g 16 ! 17 ' 27 is crucial to describe properly the low- 
energy physics of FeSe. Including only density-density 
interactions drastically suppresses the quasi-particle at 
low-energies (see Fig. [5]), rendering the system much more 
incoherent for the same value of parameters, as reflected 
in high scattering rates — ImE(iw + ) in the range 0.5 — 1.0. 
In this case also, lowering the temperature (here, down 
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FIG. 5: Difference of the total DOS when calculated using the 
full 4- index [/-matrix (black) and density-density interactions 
(red). For both calculations: U = 4.06 eV, J = 0.91 eV. 



to T = 50 K) does not reduce the scattering rates. An 
increase of the Coulomb parameters by even 50% does 
not drive the system to a Mott-insulator: at U = 6.0 eV, 
J = 1.3 eV there is still finite spectral weight at the Fermi 
level. We increased U further to U = 10 eV which fi- 
nally leads to insulating behavior. Taken together, our 
results for both rotationally invariant and density-density 
interactions reveal that there is region of the (U, J) pa- 
rameter space, where the system is in a state which is 
neither a Fermi-liquid, nor a Mott insulator. This region 
in parameter space looks similar to the selective localisa- 
tion found in Refs. mm, although we do not find one 



Mott-localised orbital, but instead three incoherent or- 
bitals. This is also consistent with recent reports within 
DFT+DMFT calculations based on exact diagonalisation 
for LaFeAsQi 7 - and FeSe.-^ 

This suggests that FeSe is just beyond the transition 
from the Fermi-liquid phase to an intermediate bad- 
metallic phase, and could be a realization of the 'spin- 
freezing' scenario . 15 ' 16 Indeed, the response to an ex- 
ternal magnetic field increases significantly, when this 
phase is entered (for instance by tuning J), whichs re- 
veals the formation of localised moments. This appears 
to be consistent with the rather large values of the resis- 
tivity measured even on single-crystal samples^ and with 
NMR^i and neutron scattering experiments^ suggestive 
of strong spin fluctuations. 



IV. CONCLUSIONS 

We have presented strong theoretical and first- 
principles evidence for enhanced electronic correlation ef- 
fects in FeSe. We have identified a Hubbard satellite in 
the spectral function, large orbital-dependent mass en- 
hancements ranging from 2 to 5, and -importantly - large 
quasi-particle damping effects that may render a Fermi 
liquid description of the normal state of FeSe invalid over 
a wide temperature range. Concerning the momentum- 
resolved spectral function, we find very good agreement 
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with recent ARPES measurments. In addition to the 
mass enhancement, we found that a correlation induced 
crystal-field splitting is important to describe size and 
position of the hole and electron pockets at the F and 
M point of the Brilloiun zone, respectively. As argued 
in Ref. H, a simple scaling of the LDA bands is not 
enough to account for the measured electronic structure. 
The present orbital-sensitive non-Fermi liquid behavior 
further indicated by their sensitivity to the Hund's cou- 
pling may serve for understanding the interplay between 
the spin-charge-orbital dynamics and superconductivity 
in iron superconductors. 
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